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Abstract 

In this paper, we study systems of nonlinear second-order variational inequalities with interconnected 
bilateral obstacles with non-local terms. They are of min-max and max-min types and related to a mul¬ 
tiple modes zero-sum switching game in the jump-diffusion model. Using systems of penalized reflected 
backward SDEs with jumps and unilateral interconnected obstacles, and their associated deterministic 
functions, we construct for each system a continuous viscosity solution which is unique in the class of 
functions with polynomial growth. 
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1 Introduction 

During the last decade optimal stochastic switching problems have attracted a lot of research activity (see 

e.g. [6, 11, 13, 12, 16, 17, 18, 19, 22, 26, 30] and the references therein) in connection with their various 

applications especially in the economic and finance spheres, such as energy, etc. Comparatively, switching 

games, of zero-sum or nonzero-sum types, have been less considered even though there are some works in 

this field including [10, 23, 24, 25, 28]. In these latter articles, the Hamilton-Jacobi-Bellmans-Isaacs PDE, 

which is of min-max or max-min type, associated with the zero-sum switching stochastic game is studied 

from the point of view of viscosity solution theory. The probabilistic version of those works is considered in 

[10, 23] where it is shown that the BSDE system associated with the zero-sum game has a solution. In [10], 

uniqueness of the solution, which is an involved question, is proved as well. The issue of existence of a value 

or a saddle-point for the game is also addressed in [10], where it is shown that the game has a saddle-point 
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when the switching costs and utilities are decoupled. This existence is deeply related to the optimal policy 
of a standard optimal switching problem. The general case still open. 

Except articles [20, 16], all the previous works deal with the case of Brownian noise. In [20], the framework 
where the noise is driven by a Levy process is studied in detail. The main motivation is that models which 
include Poisson jumps have the feature to be more realistic since they capture non-predictable events, e.g. in 
the energy market, jumps of the prices due to sudden weather changes, etc. Therefore the main objective of 
this work is the extension to the model with jumps of the paper [10], where the authors have studied systems 
of variational inequalities with interconnected lower and upper obstacles, which arise as the Hamilton-Jacobi- 
Bellman-Isaacs equation in a multiple modes switching game between two players in the framework without 
jumps. Precisely we consider the following system of non-local variational inequalities or integral-partial 
differential equations (IPDEs for short): For every pair (i, j) in the finite set of modes A 1 x A 2 , 

I min{ (v 1J - L %3 [(v H )^i) eA i xA 2 .])(t,x), max{(u y - U l3 [(v kl )^ kt i) eA i xA 2 ])(t,x), 

< -d t v ij (t,x) - Cv ij (t,x) - g ij (t,x,(v kl (t,x)) {k ' l)GA i xA 2 ,a(t,x) T D x v i3 (t,x), Iij(t,x,v ij ))}} = 0; (1-1) 

I v ij (T,x) = h ij (x) 

where, for any (t,x) £ [0, T] x R fc , 

a) C(j>{t,x) := b(t,x) T D x <j)(t,x) + \Tr[a(j T (t,x)D 2 x (f)(t,x)}+ 

+ /3(x, e)) - cj)(t,x ) - D x <j>(t,x)/3(x, e))n(de)-, 

b) Iij(t,x,<j)) = f E (c/>(t,x + /3(x,e)) -<j>(t,x))Y 3 (x,e)n(de); (1.2) 

c) L i3 [(v kl ) (k i )eA i xA2 ])(t,x) :=max(v kj - g )(f, x)) and 

k^i ~ lK 

U ij [(v kl )^ k ,i)eA 1 xA 2 ])(t, x) := min (v il +g jl )(t,x)). 

The function g (resp. Tjji) stands for the switching cost of the maximizer (resp. minimizer) when she makes 
the decision to switch from mode i to mode k (resp. mode j to mode l) while the function g IJ (resp. h lJ ) is 
the instantaneous (resp. terminal) payoff when the maximizer (resp. minimizer) chooses mode i (resp. j). 
The non-local terms which appear in (1.1) and given in a), b) above stem from the jumps of the dynamics of 
the system which is of jump-diffusion type (see (2.5) below). Finally note that the obstacles in (1.1) depend 
on the solution. 

In this paper we show that system (1.1) has a continuous solution in viscosity sense which is moreover 
unique in the class of functions which have polynomial growth. As a by-product we obtain the same conclusion 
for the max-min system (4.26). Our work should be seen as a starting point for future research in this field 
(e.g. improvement of the results, numerics, etc.). 

This paper is organized as follows: In Section 2, we fix the notations, assumptions, definitions and 
set up accurately the problem. In Section 3, we prove a comparison result between the subsolutions and 
supersolutions of system (1.1) when they have polynomial growth. As an immediate consequence, the 
solutions of (1.1) with polynomial growth is necessarily continuous and unique. In Section 4, we introduce 
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systems of integral-partial differential equations with lower (resp. upper) interconnected obstacles which are 
obtained by the penalization of the upper (resp. lower) obstacles of system (1.1). They are approximating 
schemes for (1.1) and max-min system respectively. We highlight some of their properties in making the 
connection with systems of reflected BSDEs with lower (resp. upper) obstacles. Later on we show that 
system (1.1) has a subsolution and a supersolution as well. Finally by Perron’s method we show that it has 
a unique solution. As a by product we show also that the max-min system has a unique solution. At the 
end, there is an Appendix, where we give another definition of the viscosity solution of system (1.1) which 
uses “local” maxima and minima and which is inspired by the work by Barles-Imbert in [4]. 

2 Preliminaries 

Let (0, T, {T t ) t >o, P ) be a stochastic basis such that Jo contains all P-null sets of P - , and P)+ = f) Tt+e 

£> 0 

P), t > 0. We suppose that the filtration is generated by the following two mutually independent processes: 

- a d-dimensional standard Brownian motion ( B t )t>o ; 

- a Poisson random measure N on R + x E, where E = R ( — {0} (l > 1 fixed) is equipped with its Borel 
field Be, with compensator v(dt,,de) = dtn(de), such that {lV((0,t] x A) = (N — i/)((0 ,t\ x A)} 0 <t<T is an 
Pi-martingale for all A G Be satisfying n(A) < oo. The measure n is assumed to be u-finite on (E, Be) and 
integrates (1 A \e\ 2 ) e& E- 

Let T be a fixed positive constant and let A 1 (resp. A 2 ) denote the set of switching modes for player 
1 (resp. player 2) whose cardinal is mi (resp. 7712 ). The set A 1 x A 2 will be sometimes simply denoted 
by T. For ( i,j ) € T, we set A) := A 1 — {i}, A 2 := A 2 — {j} and = F — {(i, j)}. Next, for 

V = ( y kl )(k,i)eA 1 xA 2 € R miXr ” 2 and yi G R, we denote by [y l -’,yi] the matrix obtained from y by replacing 
the element y lJ with y\. 

A function $ : (t,x) G [0, T] x R fc 1 —> $(t, x) G R is called of polynomial growth if there exist two 
non-negative real constant C and 7 such that for any ( t,x ) G [0,T] x R fc , 

mt, x) 1 << 7 ( 1 + 140 . 

Hereafter, this class of functions is denoted by n s . 

We now define the probabilistic tools and sets we need later. Let: 

(i) V be the er-algebra of Pi-predictable subsets of H x [0,T]; 

(h) H 2 := {if := ( 77 ) t<r is an Revalued, P)-progressively measurable process s.t. |MI ^2 := e (/ 0 T \Vt\ 2 dt) < 
00 } ; 

(iii) <S 2 := {£ := (^ t)t<T is an R-valued, P)-adapted RCLL process s.t. ||^||| 2 := E[sup 0<t<T |^ t | 2 ] < 00 } ; 
A 2 is the subspace of S 2 of continuous non-decreasing processes null at t = 0 ; 

(iv) , H 2 {N) := {U : fl x [0, T\ x E —> R, 'P^Hs-measurable and s.t. ||H||^ 2 ^ := E(J Q T f E \U t (e)\ 2 n(de)dt) < 
00 }. 
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The main objective of this paper is to investigate the problem of existence and uniqueness of a viscosity 
solutions v(t,x) := (V J (t, x))^j^ A i xA 2 of the following system of non-local variational inequalities (SVI in 
short) or IPDEs with upper and lower interconnected obstacles: V(i, j) £ A 1 x A 2 , 

min{(v u — L z i[v\)(t, i); max{(i) ,J ' — U l ^[v\){t,x)\ 

-d t v ij {t,x) - Cv ij (t, x) - g lJ (t, x, {v kl {t,x)) {k ^ A i xA 2 ,a{t,x) T D x v i: >{t,x),Iij(1:,x,v i:i ))}} = 0; 

(^• 1 ) 

v ij (T,x) = h ij (x), 

where for any ( i,j ) £ T, ( t,x ) £ [0,T] x R fe and (j> £ C 1 ’ 2 , L^[v\, U l i[v\, Cv l ^{t,x) and Iij(t,x,(j )) are given 
in (1.2). The functions g ik ,~gn-, /? and ' y lJ are given and will be specified more later. 

Next for S > 0, (t, x) £ [0,T] x R fc , ^ £ R fc , </> a C 1,2 -function and (i,j) £ T, let us set: 

(a) lj(t, x, 4>) = f\ e \< s {(j>(t, x + 0(x, e)) - cj)(t, x) - D x (j>{t , x)/3(x, e))n{de)\ 

(b) I$(t, x, (j)) = f\ e \> s (<f>{t,x + (3{x,e)) - <j>{t,x) - C P{x,e))n{de)\ 

(c) x , <t>) = /| e |< 5 (^(i, * + 0(x, e)) - 0 (i, x))Y J (x, e)n(de); 

(d) 0) = i|e|> 5 (^(^ x e )) - a; ))7 y (*, e)n(de); 

(e) C(j)(t,x) := b(t,x) T D x <f>(t,x) + ^Tr[crcr T (f, x)D 2 x <f>(t : a:)]. 

Note that for any 6 > 0 and (i,j) £ A 1 x A 2 , 

I{t,x,4>) = Ig(t,x,<j>) + I$(t,x,D x <t>,(j)) and Iij(t,x,(t>) = l]f{t,x,<t>) + I?/(t,x,<j>). 

Next the following assumptions will be in force throughout the rest of this paper. 

(AO) : 

(i) The function b(t,x) (resp. a(t,x)): [0,T] x R fc —► R fe (resp. R fex d ) is jointly continuous in (t, x) 
and Lipschitz continuous w.r.t. x, meaning that there exists a non-negative constant C such that for any 
(t,x,x') £ [0,T] x R fc+fc we have: 

| <r(t, x) — cr(t, x')\ + | b(t, x) — b(t, x')\ < C\x — x'\. 

Combining this property with continuity one deduces that b and er are of linear growth w.r.t. x, i.e., 

\b{t,x)\ + \a(t,x)\ < <7(1 + |a;|). 

(ii) The function j3 : R fe x E —> R fc is measurable, and such that for some real K, 

\P(x, e)| < K(1 A |e|) and |/3(x, e) — (3{x', e)\ < K\x — #'1(1 A |e|), Ve £ E and x, x' £ R fe . 

(Al) : 

For any (i, j) £ T, the function g l i : ( t,x,y,z,q ) £ [0, T] x R fc + m i xm 2 +d+i ,—^ (t,x,y, z,q) £ R verifies: 

(i) it is continuous in (t,x) uniformly w.r.t. the other variables ( y,z,q ) and for any (t,x) the mapping 
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(t, x ) i-> g t3 (t, a;, 0 , 0 , 0 ) is of polynomial growth ; 

(ii) it satisfies the standard hypothesis of Lipschitz continuity w.r.t. the variables ( y,z,q ), i.e. for any 
(t,x) G [0,T] x R fc , (yi,y 2 ) G (R miXm2 ) 2 and (^,z 2 ) G (R d ) 2 , q 1: q 2 G R, it holds 

(*, 2 /i, zi,qi) - 9 l3 (t,x,y2,Z2,q2)\ < C(|yi - y 2 \ + \zi - z 2 \ + \qi - q 2 1 ), 

where \y\ stands for the standard Euclidean norm of y in R miXm2 ; 

(iii) the mapping q g l3 (t,x,y,z,q) is non-decreasing, for all fixed ( t,x,y,z ) G [0, T] x R fc + m i xm 2 +d_ □ 

Next for any (i,j) G T, the function 7 * J : R fc x£->R verifies for some constant C : 

(a) | 7 ,J (a;, e) — 7 u (x / , e)| < C\x — x'|(l A |e|), x, x' G R fe and eGfi; ^ 

(b) 0 < 7 * J (x, e) < (7(1 A |e|), x G R fc and e G E. 

Finally let us define functions on [0,X 1 ] x R fc + m i xm 2 -W x £^(E,Be,ti), as follows: 

f 3 (t,x,y,z,u ) := g lJ (t,x,y, z, J E u{e)j lJ (x,e)n(de)). (2.3) 

(A2): Monotonicity: For any (i,j) G T and any ( k,l ) ^ (j,j), the mapping (t,x,y, z,u) is non¬ 
decreasing. □ 

(A3): The non free loop property: The switching costs g. k and gji are non-negative, jointly continuous in 
(t,x), belong to II S and satisfy the following condition: For any loop in A 1 x A 2 , i.e., any sequence of 

pairs (ii,ji), (inJn) of A 1 x A 2 such that (i N ,j N ) = card{(n, ji),..., (inJn)} = N - 1 and 

V q= 1,..., A — 1, either i q+1 = i q or j q+1 = j q , we have: 

V(t, x) G [0, T] x R fc , Y, 'Pi'j'fa*)? 0 ’ (2- 4 ) 

q=l,N—l 

where for any q = 1,..., N — 1, g> iqjq (t, x) = -9 iqiq+i (t, x)\ iq jti q+1 + g Jqjq+1 ( t , x)l jq ? jq+1 . □ 

(A4): The functions hA : R fc —> R are continuous w.r.t. x, belong to class II S and satisfy: 

V(«,j) G F, x G R fc , max(h kj (x) -g. k (T,x)) < h l3 {x) < min {h l \x) + g jl (T,x)). □ 

To begin with let us point out that the non-local terms X(t,x,<j>) and Iy(f,x, </>) introduced previously are 
well defined under Assumptions (A0) since for any function <f> of class C 1, , by the mean value theorem, we 
have 

\(j>{t,x + (i{x,e)) - <j)(t,x) - D x c/)(t,x)P{e,x)\ < C' f ( > 1 j|/3(x,e)| 2 < C^J( 1 A |e|) 2 , 

and 

|7 i3 (x,e)(0(t,x+/3(x,e)) -<j>(t,x)) \ < cj 2 J\/3(x,e)Y 3 (x,e)\ < C t (2 J( 1 A |e| 2 ) 
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where Cjf 1 ,) and Cj 2 J are bounded constants. They are the bounds of the first and second derivatives of 
y <j)(t, y) in B(x, Kp) where Kp is a bound of the function f3. □ 

Let us consider now the following SDE of jump-diffusion type ((t,x) G [0,T] x R fc ): 

x t,x = x + f t s b(r, X( ,x )dr + // a(r, Xp x )dW r + // f E /3(X*’f, e)N(dr, de), s € [t, T]. (2.5) 

The existence and uniqueness of the solution X t ’ x := (V* ,x )se[t,Tl follows from [15]. 

We now precise the definition of the viscosity solution of system (2.1). First, for a locally bounded function 
u: (t,x) G [0,T] x R fc i-)- u(t,x) G R, we define its lower semi-continuous (lsc for short) envelope u*, and 
upper semi-continuous (use for short) envelope u* as following: 

u*(t,x) = lim u(t',x'), u*(t,x) = lim u(t',x'). 

t’<T (t',x')->(t,x), t'<T 

Definition 2.1. A function u(t,x) = (u l i(t, x))^j' )e A 1 xA 2 ■ [0, X 1 ] x R fc — > R m i xm 2 such that u l i is lsc 
(resp. use) and belongs to T\ g , is said to be a viscosity supersolution (resp. subsolution) of (2.1) if for any 
(i,j) £ A 1 x A 2 : 

- u l i(T,x o) > (resp. <) h^{x o),Vxo € R fc . 

- for any (to,Xo) G (0,T) x R fc and any test function f> G C 1 , 2 ([0,T] x R fc ) such that (to,Xo) is a global 
minimum (resp. maximum) point of u lJ — <j) and u lJ (to,Xo) = 4>{to,xo) , 

min{(u u - L lJ [u\)(t 0 ,x 0 ); max{(tt ’ 3 - U l3 [u])(t 0 ,x 0 )\ 

< -d t (j>(to,xo) - b(t 0 ,x 0 ) T D x (j)(t 0 ,x 0 ) - \Tr[aa T (t 0 ,x 0 )D 2 xx (t>{to,x 0 )] - I(to,xo,<j>) 

-g v {to,xo, (u kl {t 0 ,x 0 ))(k,i)eAixA 2 ,v T (to, x 0 )D x <j>(to,xo),Iij(to,x 0 , <£))}} > 0 (resp. < 0). 

A function u = (u lJ (t, ^))(i,j)e J 4 1 xA 2 of n g is called a viscosity solution of (2.1) if (uf (t, ®))(»,j)eAi xA 2 (resp. 
(u*j(t, x))(.j j)ca 1 xA 2 ) is a viscosity supersolution (resp. subsolution) of (2.1). 

Remark 2.1. By taking gji = +oo (resp. g = +oo) for any j,l G ^4 2 (resp. i,k G A 1 ) we obtain the 
definition of a viscosity solution of the system of variational inequalities with interconnected lower (resp. 
upper) obstacles. □ 

3 Uniqueness of the viscosity solution of the non-local SVI (2.1) 

In this section we will show the uniqueness of the viscosity solution of (2.1) as a corollary of a comparison 
result. In the same way with [9], Lemma 4.1, we can prove the following lemma. 
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Lemma 3.1. Let (u zj )^j)^a 1 xA 2 (resp. (w ZJ )(i,j)eA 1 xA 2 ) be an use subsolution (resp. Isc supersolution) of 
(2.1) which belongs to II g . For (t,x) £ [0,T] x R k , let r(f, a:) be the following set: 

T(t,x) := {( i,j ) £ A 1 x A 2 ,u 13 (t,x) — w 13 (t,x) = max (u kl {t,x) — w kl (t.,x))}. 

(k,l)eA 1 xA 2 

Then there exists (io,jo) £ r (t,x) such that 

u lo3o (t,x) > L lo3 ° [«] (t, x) and ro' o7 ° (t,x) < U lo3o [w\(t,x). (3.1) 

We now give the comparison result. 

Theorem 3.1. Let u = {u ZJ )^j-\ e A 1 xA 2 (resp. w = (w ZJ )/i^ e A 1 xA 2 ) be an use subsolution (resp. Isc 
supersolution) of (2.1) which belongs to II g . Then it holds that: 

V (i,j) £ F and (t,x) £ [0,T] x R fe , u l3 (t,x) < w l3 (t,x). 

Proof. Let C and p be positive constants, which exist thanks to the polynomial growth of u and w, such 
that for any (i,j) £ T, 

\u lj (t,x)\ + \w lj (t,x)\ < (7(1 + \x\ p ). 

There exists a positive constant Ao such that for any A > Ao and 9 > 0, u e := (u 13 (t, x)—6e xt ( l+|a;| 2p+2 )( i> j) g r 
(resp. w e = (w l3 (t,x) — 0e xt ( 1 + |a;| 2p+2 )(i, g )er) is a subsolution (resp. supersolution) of (2.1) (see [20], 
pp.1634). Therefore it is enough to show that u e < w e and to take the limit as 9 —> 0 to obtain the desired 
result. Finally with the previous statement w.l.o.g one can assume that there exists a real constant R > 0 
such that for any |x| > R , w l3 (t,x) > 0 (resp. u Z3 (t,x) < 0) for any (■ i,j ) and t £ [0,T]. 

The proof now will be divided into two steps. 

Step 1: Let C'ij be the Lipschitz constant of g 13 w.r.t. y. We first assume that there exists a constant 
Ai > m 1 m 2 -J 2 (i,j)erC'ij such that for any (i,j) £ T, 

9 l3 (t,x, [y 11 , Cl], z,q) - g lJ {t,x, [f 3 , C 2 ], z, q) < -A(Ci - C 2 ), (3.2) 

for any £1 > (2 in R and (t, x, y. z , q) in their respective spaces. 

We proceed by contradiction. Assume there exists (i, x) £ [0,T] x R fc such that: 

max ( u 13 — w l3 )(t,x) > 0. (3.3) 

(*d)€ r 

Therefore there exists (t*,x*) £ [0, T[xB(0, R) (B(0,R) is the open ball in R fc centered in 0 and of radius 
R and w.l.o.g. we assume t* > 0) such that: 

max max Uu Z3 — w l3 )(t, x)} = max max Uu 13 — w 13 )(t, x)} = max {{u 13 — w l3 )(t*, a;*)} > 0 . 

Me[o,r]xR‘(»,i)er“ (t,i)e[o,r]xB(o,fi) (*,j)er (»,j)e 

(3.4) 




Next let (io,jo) € T{t*,x*) that satisfies (3.1). For e > 0, p > 0, let be the function defined as follows: 


K°t°{t, x ,y) := (u lo3o (t,x) -w lo3o {t,y)) - 


\x - y | 2 


- |t-t *| 2 -p\x--. 


Let (to,xo,yo) be such that 

^°p°(*o,zo,yo) = max 

(t,x,i/)€[0,T]xB(0,K) 2 (t,i,t()6[0,T|xR*+ i 

The second equality is valid since when \x\ > R (resp. \y\ > R), u lo3o (t,x) < 0 (resp. w lo3o (t,y) > 0). On 
the other hand (to,Xo,yo) depends actually on e and p which we omit for sake of simplicity. Next as usual 


we have 


lim(t 0 ,x 0 ,yo) = (t*,x*,x*), lim —-— = 0 , 

e —>0 e ->0 £ 


lim {u’ o3 ° (t 0 , x 0 ), w lo3 ° (t 0 ,y 0 )) = ( u l ° 3 ° (t*,x*),w l ° 3 ° (t*,x*)). 

e->0 


Therefore for e small enough it holds 


u lo3o (t 0 ,x 0 ) > max (u kj0 (to, x 0 ) — g. (t 0 ,x 0 )), 
keA) — l ° K 


w l ° 3o {to,yo) < min {vf° l {t 0 ,yo) + 9 jo i(to,yo))- 


Once more for e small enough, we are able to apply Jensen-Ishii’s Lemma for non-local operators established 
by Barles and Imbert ([3], pp. 583) (one can see also [5], Lemma 4.1, pp. 64) with u lo3 °, w l ° 3 ° and 
(j>(t, x, y) := + \t — t*\ 2 + p\x — x*\ 4 at (t 0 , x 0 , yo). For any S € (0,1) there exist p e u , q e u in R, p e w , q. € w in 

R fe and two symmetric non-negative matrices of R fcxfe such that: 

(i) 

P € u~Pw = dt<t>{to, x o, yo), Qu = dx(t>{to,Xo, yo), q e w = -d y (j>(to,x 0 , yo) ( 3 - 9 ) 


K 0 \ < (DlMto,x 0 ) 0\ + 4 4 -l 

0 -M ?,) ~ 0 0 / e -I k I k 


where x) := p\x — x *| 4 + \t — t* | 2 ; 


(ii) -p e u -b(to,xo) T ql - -Tr(cr(to,xo) M*cr(to, x 0 )) - g lojo (to, x 0 , (u lJ (to, <&(>))(*,A 2 , °(to,xo)q e u , 

!lfj 0 {to, x 0 ,(Kto,.,yo)) + lff jo {to, x 0 , u 1030 )) - l]{t 0 , x 0 , <t>(t 0 ,., y 0 )) - I 2 {to, x 0 , q e u ,u 1030 ) < 0; 

( (3-11) 

(iii) -Pw - b {to,yo) T ql - ^Tr{a{t 0 : y 0 ) T Mi u a{t 0 ,y 0 )) - g 1030 (t 0 , yo, {w lJ (t 0 , Vo))(i,j)eA^xA^, a {to, Vo) T q e w , 


+ lff in {t 0 ,x 0 ,w 1030 )) - Ig{to,yo,-4>{to,xo, ■)) - Is{to,yo,q e w ,w lo:>0 ) > 0. 


(3.12) 
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Next we are going to provide estimates for the non-local terms. By the same argument as in [20] pp.1645, 
we have: 

Is(t 0 ,x 0 ,q e u ,u iojo ) - I$(t 0 ,y 0 ,q e w ,w iojo ) < C lx °~ e Vo 1 Jj e |> g (l A \e\) 2 n{de) + I^(t 0 ,x 0 , D x ip p (t 0 ,x 0 ),ip p ) 

< C 1 ^-^ 1 + Is( t o,x 0 ,D x ip p (t 0 ,x 0 ),ipp). 

(3.13) 

On the other hand 

D 2 XX W, x, y) = 2e~ 1 I k + D 2 xx ip p (t, x, y), D 2 yy W, x, y) = 2e~ 1 I k 
and by Taylor’s expansion we have 

W, x + fi(x, e), y) — W, x) - D x <j>(t, x, y)/3(x, e) = f*(l - t)fi(x, e) T D 2 XX W, x + tfi(x, e), y)fi(x, e)dt. 

It implies that 

x 0 , Wo,., yo)) = f je |< 5 n(de) f 0 (1 - t)fi(x 0 , e) T D 2 x Wo, x 0 + t/3(x 0 , e),y 0 )fi(x 0 , e)dt 

and similarly 

Is(to, Vo, -<p(to,xo, •)) = - /| e |< 5 n{de) f 0 (1 - t)/3(y 0 , e) T D 2 y Wo, x 0 , yo + tp{y 0 , e))fi{y 0l e)dt. 
Consequently it holds that 

lim^o l}{to,x 0 ,Wo, -,yo)) = lim5-> 0 /j(to,yo, ~Wo> x o, .)) = 0. (3.14) 

Next by the definition of (to,Xo,yo), we have 

u iojo {t 0 , x 0 + /3(x 0 , e)) - u io3o (t 0 , x 0 ) 

< w z °’ 3o (t 0 ,y 0 +/3(y 0 ,e)) - w lo3o (t 0 , y 0 ) +ip p (t 0 ,x 0 + P(x 0 ,e)) - ip p (t 0 ,x 0 ) 
+W{\P(x 0 ,e) — /3(yo,e)\ 2 - 2(x 0 - yo)(P(x o,e) - p(y 0 ,e))}. 

Since 7* M0 is nonnegative, and by the assumptions on /3 (see (A0)-(i)), for any S > 0, 

^ojo^O’Xo, u l0j0 ) - I?’ S jo (to,yo,w* o:i0 ) < I? o S jo (t 0 ,x 0 ,ip p ) + O(e~ 1 \xo - y 0 \ 2 ), 
and it is easy to check that 

l4bo(*o,zo, Vv)l < C Pf\z\>s 1 A |e| 2 n(de). 

On the other hand, since <f> is a C 2 -function then, using once more Taylor’s expansion to obtain: 

4ojo (to, x o, Wo, ; yo)) = f\ e \< 5 {Wo, X 0 + (3{x 0 , e),y 0 ) - Wo, x 0 , yo)}7 loJ °( x o, e)n(de) 


(3.15) 


(3.16) 


= L l<s n(de) f 0 dt{ 1 - t)/3(x 0 , e) T [§{x 0 + t/3(x 0 , e) - y 0 } + D x ip p (t 0 , x 0 + t/3(x 0 , e))\W° 0 [x 0 ,e) 


xp-yo 

e J |e|< 


/|e|<<5 P( x o, e) T 7 ioJ °( xq, e)n(de) + § /, e ,< 5 |/3(xo,e)| 2 7 io ^(x 0 ,e)n(de 


+ i | e |<5 n( de )P( x o, e) T DWWo, xo + tp(x 0 , e)) 7 l ° 30 (x 0 , e) 


(3.17) 
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and 

l!; S jo (to,yo,-<t>(to,x 0 ,.)) = ~ f\ e \<s{<t>(to,xo,yo + /%o,e)) - <f>(t 0 , x 0 , yo)}Y° 3 °(yo, e)n(de) 


which implies that 


f\e\<s P(yo, e) T 7 iw '°(j/ 0 , e)n(de) - § f M<g |/?(x 0 , e)| 2 7 iw ' o (x 0 , e)n(de) 


|e|<(5 


(3.18) 


lim^ 0 /M o (to,x 0 , <£(*(),-,yo)) = lim^o jM o (t 0 , 2 /o,-?!>(io,xo,.)) = 0. (3.19) 

Making now the difference between (3.11) and (3.12) yields 

~(Pu ~Pw) ~ [b(t 0 ,x 0 ) T q e u -b(t 0 ,yo) T qtn\ ~ \{ Tr W(to, x 0 ) T M*a(t 0 , x 0 ) - a(t 0 ,y 0 ) T M^a(t 0 ,y 0 )}} 

~[g iojo (t 0 , x 0 , (u h (t 0 ,xo))(i,j)eAixA2,(7(to,xo) T ql, lj’ s jo (t 0 , x 0 , <t>(to,., yo)) + lff jo (to, x 0 , u io3 °)) 

-g ioj °(to,yo, (w ij (t 0 ,yo))(ij)eAixA2,cr(to, yo) T q e w , ll’j 0 (to,yo, ~4>(to, Xo, ■)) + I?’j 0 (to,yo, w iojo ))} 

- I l s (to,xo,(t)(to , -,yo)) + I 1,s (to,yo, —<t>(to, xo,.)) - I 2 ’ 5 (t 0 ,xo,qu,u iojo ) + I 2 ’ s (t 0 ,yo,qw, w io3 °) < 0. 

(3.20) 

As q G R K > g 13 (t,x,y, z,q) is non-decreasing and Lipscliitz then, by linearization procedure of Lipschitz 
functions, there exists a bounded non-negative quantity H (which depends on e, 6, etc.) such that 

- [g loJO (to,x 0 , (u lJ (to, x 0 ))(ij) eA ^xA^cr(to,Xo) T q e u ,llf jo (to, x 0 ,(/)(to, .,y 0 )) + I?f jo (t 0 ,x 0 , u l ° 30 )) 

- g l ° 3 °(to,yo, (w' 3 (to,yo))(i,j)eAixA*,(r(to,yo) T qli),llfj 0 (to,yo,-4>(to,xo ,.)) + tffj 0 (to,yo,w l ° 3 °))\ 

= -[ g 1030 (to, Xo, (u 13 (to, x 0 ))(ij)eAi xA^,(x(to,xo) T q e u, llfj Q (to, xo, 4>(to, ■, Vo)) + lff jo (t 0 ,x 0 , u l ° 30 )) 

- g ,OJO (to, yo, (w ,J (to, yo))(i,j)GA i X A^,cr(to,yo) T q e w ,ll’j 0 (to, x 0 , <j>(t 0 ,., y 0 )) + /f 0 j 0 (t 0 , x 0 , u Wo ))] 

- s x Vlfj 0 (to, xo,<t>(to, ; yo)) + lff jo (t- 0 , Xo, u 1030 ) - lh s jo (to, yo, -<j>(t 0 ,x 0 , ■)) - lff jo (t 0 , yo, w 1030 )} 

On the other hand, once more by a linearization procedure and taking into account of (3.2), the first term 
of the right-hand side of the previous equality verifies: 

[g lo3o (t 0 ,x o, (u 13 (to,xo))(ij) eA 'xA 2 ,(T(to,xo) T q e u , llf jo (to,x 0 ,(l>(to, -,yo)) + lff jo (to,x 0 ,u l ° 30 )) 

- g l03 ° (t'O, yo, (W 13 (to, yo))(i,j)eAi xA*,<r(to,yo) T q e w ,llf jo (to, X 0 , </>(t 0 , ., yo)) + lff jo (to, Xo, ■U W0 ))] 

< -\(u lo3o (to,x 0 ) - W lo3o (t 0 ,yo)) + ^2 Ct 0 , X0 ,y 0 ,s(u lJ (to,x 0 ) - W 13 (t 0 ,y 0 )) 

(*>j)^(*odo) 

+ Vt 0 0 J ,xo,y 0 ,s( (J (to,xo) T q e ~ cr(to,yo) T q e w ) + SU P W o3 °(to,xo,y,z,q) - g lo3 °(to,yo,y,z,q)\ 

V,z,q 

where ^ XQ yo s is bounded non-negative quantity (positiveness stems from (A2)) by Ci 0 j 0 , and rj l t J Q xq yo § is 
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a bounded quantity by the Lipscliitz constant of g l °i° w.r.t. z. Therefore from (3.20) we deduce: 
X(u iojo (to,x 0 ) - w iojo (t 0 ,yo )) < 

( Pu-Pw ) + [b(to,xo) T q e u - b{to,yo) T q e w \ + ±{Tr[a(t 0 ,x 0 ) T M^a(t 0 ,x 0 ) - a(t 0 ,y 0 ) T M^a(t 0 ,y 0 )]} 

+ s -[ J * 1 0 io( t o> a: o,^(io,..,yo)) ~ llfj 0 {to,yo,~(l>(to,xo,.)) + lff jo (to,x 0 ,u iojo ) - lf’ s jo (t 0 ,yo,w io:io )] 

+ I x ' 5 (to,xo,(j){to,.,yo)) - I'^ito^o^dito^o,.)) + I 2 ’ 5 (to,x 0 ,qu,u iojo ) - I 2 ’ S (t 0 ,yo 1 qu„w iojo ) 

+ E(i,j)^(* 0 ,j 0 ) Cio^o.OT^^o^o) - W ij {t 0 , y 0 )) + ??t°f° 0>OT (cr(to,a:o) T g' ! - cr(<o,yo) T ^) 

+ su Py,z,« |5' ioio (^o,a;o,y,2,g) - 9 iojo {t 0 ,yo, y,z,q)\. 

(3.21) 

Next recall (3.13), (3.14),(3.15)-(ii) and (3.19), then take the limit superior as 6 —> 0 then the limit superior 
as p —» 0 to obtain: 

X(u iojo (to,x 0 ) - w iojo (t 0 ,yo )) < 

(p« -P^) + [^(^o,a;o) T g e u - t»(t 0 , 2 /o) t 9^] + ^{Tr[a(to,xo) T M*a{to,x 0 ) - a(t 0 ,y 0 ) T M^a(t 0 ,y 0 )]} 

+ E(i,j)#(i 0 ,j 0 ) Cy(w ij (io,a:o) - ^(t 0 ,2/o)) + + C iM o(cr(f 0 ,2;o) V - o-(t 0 ,2/ 0 ) T <7^) 

+ sn Py,z, q \ 9 iojo (to,x 0 ,y,z,q) - g iojo (t 0 ,y 0 ,y,z,q) \ + |0(e _1 |a;o - 2/o| 2 )| + C'e _1 |a:o - Vo\ 2 - 

(3.22) 

Finally by the continuity of g l °i° (see (Al)-(i)), using the properties satisfied by p* q e u , q^, M* and M£, 
sending e to 0 and taking into account of (3.5)-(3.6) to obtain (in a classical way) that: 

X{u iojo (t*,x*) - w iojo (t*,x*)) < Cij(u ij {t*,x*) - w ij (t*,x*))+ 

which is contradictory by the definitions of A and (t*,x*). Thus for any (i,j) £ A 1 x A 2 , id- 7 < w l T Note 
that we have used the fact that id- 7 (resp w lJ ) is use (resp. lsc) when we take the limit as e —» 0 to deduce 
the last inequality. 

Step 2: The general case. If (id-^^jjg^i xA 2 (resp. (ud- 7 )^^^! xA2 ) is a subsolution (resp. supersolution) of 
(2.1) in the class II S , then for any A £ R, {e~ xt u i: >)^ i ^ eA i xA 2 . (resp. (e~ Xt w i: ’)( i j' )eA i xA 2 ) is a subsolution 
(resp. supersolution) of a system of type (2.1) but associated with 
{(e _A V y (*, x, {e Xt u kl )( k j ) eA i xA 2 , e xt a(t, x) T D x u i: > (t, x), I y (i, a;, e xt u lJ )) - Aid- 7 (f, a:) 

( e_At £j/ s ( i ’ a: ))bfee J 4 1 > ( e ~ Xt 9jl(t,x))jj e A 2 , {e~ XT h %3 {x))(ij) eA i xA 2 )}. But in choosing appropriately the con¬ 
stant A we get that the functions g^(t, x, y, 2 d- 7 ,q) := e~ Xt g i: >(t, x, e xt y , e xt z l \ q) — Ay 1 - 7 verify (3.2). Therefore 
by the result of Step 1 we have for any (i, j) £ A 1 x A 2 , e~ Xt u'^ < e~ xt w 1 ^ , whence the desired result. □ 

As a by-product we have: 

Corollary 3.1. The system of variational inequalities with bilateral obstacles (2.1) has at most one viscosity 
solution in the class II 9 which is moreover necessarily continuous. 


Finally a remark in the case when we have only lower or upper interconnected obstacles. 
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Remark 3.1. If we assume that gjk = +oo (resp. = Too,) for any j,k £ A 2 (res. i,l £ A 1 ), then 
system (2.1) turns into of type with one lower (resp. upper) interconnected obstacles. As the non free loop 
property is satisfied by (resp. ( 9jk)j,keA 2 ) then the system of variational inequalities with lower 

(resp. upper) interconnected obstacles has at most one viscosity solution in the class II g which is moreover 
necessarily continuous. □ 


4 Approximating schemes and BSDEs 

For n, m > 0, let (Y , ^ ,n,m , U 1 h n ’ m )^ i ^ GA i xA 2 be the solution of the following system of BSDEs: 


(Z S 2 X U 2 X H 2 (A); 

fij.n.m (~ -yt.x ( \rkl,n,m\ ryij.n.m TTii,n.m\A „ i 

dy s = -J P,A S ’ A Y s )(k,i)eA^xA^,A s J ’ ’ ,U S J ’ )as+ 
Z(d ,n ’ rn dB s + f E U i J’ n ’ m (e)N(ds, de), s < T; 

Y h),n,m = h ii( X %*), 


(4.1) 


where 




= g ij (s,X t s ’*,(y k %j )eAlxA 2 ,z,f E u( e ) ^(X^,e)n(de)) 

+ n(y lJ - max{j/ fcj - g (s, X\’ x )})^ - m.{y ij 
keA} 


f E u(e)Y J (X*’ x ,e)n(de)) 

- pi }]{y d +9ji(s,Xl’ x )}) + . 
lGA j 


Let us notice that under Assumption (Al), the solution (Y l h n < m , Z‘^ n ' m ,U‘ l ^ n ' m )^ i ^ eA i x a 2 of (4.1) ex¬ 
ists and is unique (see e.g. [29] or [3]). Next, let us focus on the properties of the matrix of processes 
A )(;,.der- 


Proposition 4.1. For any ( i,j ) £ A 1 x A 2 and n,m> 0 we have: 

(*) 

P - a.s ., Y ij ’ n ’ m < Y ij ’ n+ i’ m and Y ij ’ n ’ m+1 < Y ij ’ n ’ m . (4.2) 

(ii) There exists a deterministic continuous function v l h n,m € n g s.t., for any (t,x) € [0,T] x R fc , 


Yfi’ n ’ m = V i i’ n ’ m (s,X t s ’ x ), 


s £ [t, T], 


(4.3) 


(Hi) for any (t,x) £ [0,T] x R k , 

v ij ’ n ’ m (t,x) < v ij ’ n+1 ’ m (t,x ) and v ij ’ n ’ m+1 (t, x) < v ij ’ n ' m (t,x). 


(4.4) 


Proof, (i) Let n and m be fixed. We are going to use a result by X.Zhu ([31], Theorem 3.1) related to 
comparison of solutions of multi-dimensional BSDEs with jumps. It is enough to show that for any t £ 

[0,T},( yij ) {iJ)er , (%) ( i,j) e r e R miXm2 , {zij)(i,j)er, (zij)(i,j)e r e ( Rp‘X">= and (« !3 ) (i , j)er , (%)(*, tier € 
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(. £ft(E,13E,n ) miXm2 , there exists a constant C such that: 

-4 yij{f lJ ’ n+1 ' m ( s ^ X s X ’(yki +yki)(k,i)er,z i j,u ij ) - / u ’"’ m (s,X*’ x ,(y w )( fci/)e r,^j,Wy)} 

(*,j)er 

— 2 ^ ' 1{j/O'< 0} l^ij — Zjj I + C 'y ' {.Vij) 

(i,j)e r (yjgA'xA 2 


+ 2 / l{y«>0}l(yu +Mij(e) -u*j(e)) | 2 n(de) 

x- Vr JE 

+ 2 ^ ' / l{j/*j<o}[|(yij + Uij(e) — Uij(e)) \ — \y i j\ — 2 Uij{uij(e) — u-ij(e))]n(de) 

t ■ ■ \ ^-n £7 


(hi)er 

where for any itR, i + =iV 0 and x~ = (—x) V 0. But the above inequality follows from the facts that 

for any (i.j) £ T: 

^ j?ij,n,m ^ j?ij,n-\-l,m . 

b) For any (■ 9ki)(k,l)er £ R miXm2 such that 9 k i > 0 if (k,l) £ anc l g.j = o ; 

( y kl + 0 w )(k,j)er , Zij , Uij ) > r’ n ’ m (s,Xl’*, (y fc 0(JM)er, *«,««)• 

c) f lJ depends only on Zij,Uij and not on the other components z k i,u k i, (k , l) ^ ( i,j ). 

d) 

- 1//,,(,s-. .V.:-' . (yhi)(k,i)er,Zij,Uij) ~ f ij ’ n ' m (s,Xi’*, (y fcI )(fc,i)er, *«,*«)) 

< C(^) 2 + 2 / l{j/y>o}l(l/y + «y(e) - %(e)) _ | 2 n(de) 

JE 

+ 2 / l{yy<0}[|(2/ii + 'M e ) -%-(e))-| 2 - | y r| 2 - 2y ij (u ij (e) - u xj (e))]n(de) 

JE 

Indeed a), b) and c) are easy to check. We just need to prove d). In the case when ty, > 0, this inequality is 
obvious. Next let us focus on the case when y l: j < 0. First note that by a linearization procedure we have: 

- 4 y^(f ij ’ n ’ m (s,Xl’*, (y fcI ) (M)e r, *«,««) - f ij ’ n ’ m (s,Xl’*, (|7 w )( fclI ) 6 r,^,«tf)) 

= -4y^(fl y (s,X s t > :r ,(y fcJ ) (fc>0er ,^ i , f 7* J (X s t ’ a: ,e)u li (e)n(de)) 

J E 

/ 7 ,J (^ ,a: ,e)u*j(e)n(de))) 

JE 

= ^Vij x Si x / 7 y (X*’ x ,e)(uy(e) - t%(e))n(de) 

JE 

where Si is a non-negative quantity (since y' J is nondecreasing in g) and bounded by the Lipschitz constant 
of w.r.t. g and which depends on the other variables. Now 

4 y~j x Si x f E Y J (X*’ x ,e)(u ij (e) - wj(e))n(de) 

Y J (X t s ’ x ,e)(u ij ( e ) - Uij(e))n(de) 


= 4 y t3 x =i x 


* Uij Hij ^ Vi] 


+ 4 y t , x Si x 


7 u (^s ,e)(uij(e) - Uij(e))n(de) 


< 4y y - x Si x 


Hij ^ yij 

'Y t3 (X t a ’ x ,e)(u ij (e) - «ij(e))n(de). 


Hij Uij < yij 
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But for an appropriate constant C, 

Y J (X t s ’ x ,e)(u ij (e) - < C{^ 3 {Xf x ,e)) 2 + 2{u ij (e) - ?%(e)) 2 

then 

4 y~j x Si x f E 7 ij (Xp x ,e)(uij(e) - u xj (e))n(de) 

< C (Vij) 2+ 2 S Uij -u ii <- ViJ 1 {v i ’<o}[\(yij + u ij ( e ) -Uij{e ))~| 2 - |t/^| 2 -2y ij (u ij (e) - u*j(e))]n(de) 

< C(y~j) 2 + 2f E t{ y ij <0 }[\(yij + Uij{e) -Uij{e))~\ 2 - |%| 2 - 2j/ y - («y(e) - %(e))]n(de) 

which is the claim. Thus P - a.s., Y ij ’ n ’ m < Y ij ’ n+1 ’ m . In the same way we have also P - a.s., Y ij ’ n ’ m+1 < 

'y r ^ji n ,m 

The second claim is just the representation of solutions of standard BSDEs with jumps by deterministic 
functions in the Markovian framework (see [3]). The inequalities of (4.4) are obtained by taking s = t in 
(4.2) in view of the representation (4.3) of Y l3 ’ n ’ m by v l3,n ’ m and X t,x . □ 

Remark 4.1. Forany(i,j) ^Y,v l3 ’ n,m is the unique viscosity solution in A g of the following integral-partial 
differential equation: 

—dtv l3,n ’ m (t,x) — Cv l3 ' n,m {t,x) — n(v l3,n,rn (t,x) — rna x(v k3 ' n,m (t,x) — g (t,x))~ 

keA\ 

+m(v l3,n ’ m (t,x) — min (v ll,n ’ m (t,x) + cjji(t,x)) + 

leA 2 

a;, (^’->^(i,^)) (fe;i)e r,a;),(i,a:,^->-)) = 0; (4 ' 5) 

v ij >’ m (T» = h i3 (x). 

For more details one can see [3]. 

We now suggest two approximation schemes obtained from the sequence {{Y t 3 ,n ’ m )(i,j)er)n,m of the 
solution of system (4.1). The first scheme is a sequence of decreasing reflected BSDEs with interconnected 
lower obstacles defined as: V(i,j) £ A 1 x A 2 , 

( (Y i3 ’ m , Z i3 ’ m , U i3 ’ m , K i3 ’ m ) e «S 2 x n 2 x n 2 (N) x A 2 -, 


Y*i' m = h i3 (X t f x ) + f s J f i3 ’ m {r,Xf x ,{Y^ m ) ( ^ l)er ,Zi 3 ' m ,Ui. 3 ’ m )dr - fj Z?' m dB, 
- Is Se Uf. 3 ’ m (e)N(dr, de) + K%’ m - K i3 ’ m , s < T; 


Y, ij ’ m > m ax{F^’ m - g (s, Xf x )}, s < T; 

k£A\ ~ lK 

fo(Yj 3,m - ma x{Y s k3 ’ m - gjs, Xf x )})dK}= 0 , 


(4.6) 




15 


where for any (i,j) £ T, to > 0 and ( s,y,z,u ) (u £ C^E, Be,™)), 

f ij ’ m (s, X*’*, V, 2, u) := g ij ’+’ m (s, Xl' x , {y kl )(k, l)eA i xA 2 ,z, f E u(e)'Y ij (X*’ x , e)n(de)) 

:=g»(s,X^,(y kl ) ( k,i )eAlxA 2 ,z,f E u(e) 7 i3 (X^, e)n(de)) 

-m{y ij - min(y il +g jl (s,X t s ’ x )))+. 
1&A 2 : 


(4.7) 


The following result is related to existence and uniqueness of the solution of (4.6) and some of its properties. 
Proposition 4.2. i) For any fixed m > 0, the solution (Y* 3,m , Z Z3,m , C/*- 7 ’ 7 ", K l3,m )^^) er of the system (4-6) 
exists and is unique. Moreover for any ( i,j) and m > 0, we have: 


lim E[ sup | Yf j ’ n ’ m - Y s ij ’ m | 2 ] -> 0 and P - a.s., Y ij ’ m > Y ij ’ m+1 . (4.8) 

n->oo t<s<T 

(ii) There exists a unique mi x ni 2 ~uplet of deterministic continuous functions (u kl,m )(k,l)eA 1 xA 2 m n g such 
that, for every t <T, 

Y s ij ’ m = u ij ’ m (s, Xp x ), s £ [i, T\. (4.9) 

Moreover, for any ( i,j ) £ T and ( t,x ) £ [0,T] x R fe , u l3,m (t,x) > u l3 ’ m+1 (t,x). Finally, (u Z3,m )^jj eA i xA 2 
is a unique viscosity solution in the class II g of the following system of variational inequalities with inter¬ 
connected obstacles: V(i,j) £ T, 

?nin{u v,m (t, x) — max(u kj,m (t,x) — g (t,x)); —d t u l3,rn (t,x) — Cu l3,m (t, x) 
k£A] ~ lk 

' -g^+’ m (t,x,(u kl ’ m (t,x))^ l)er ,a(t,x) T D x u^ m (t,x),I ij (t,x,u^ m )} = 0; ( 4 - 10 ) 

u ij ’ m {T,x) = h ij (x). 

Proof, (i) It is enough to consider the case m = 0, since for any (i,j) £ A 1 x A 2 , the function 

(s,x, (y kl )(ki)eA^xA 2 ) -t -m{y 13 - min (y zl +g jl {s,x))) + 

IdA 2 

has the same properties as f 13 displayed in (Al) and (A2). First and w.l.o.g we may assume that f 13 is 
non-decreasing w.r.t. y kl , for any (k,l) £ T, since thanks to assumption (A2), it is enough to multiply the 
solution by e At , where A is appropriately chosen in order to fall in this latter case, since f 13 is Lipschitz w.r.t. 
the component y 13 . Now let FI (resp. F) be the following functions: 

H(x) = \h Z3 {x)\ and F(t,x,y,z,u) = ^ \f 3 (t,x,yl mi ,m 2 , z,u)\ 

r r 

where ( y,z,u ) £ R 1+d x C^{E,BE,n) and I mi ,m 2 is the matrix of m i (resp. m 2 ) rows (resp. columns) 
with entries equal to 1. Let (Y, Z, U) be the solution of the following one-dinrensional BSDE with jumps 
associated with (F(s, ,y, z,u), H(X E x )). Next let n be fixed and let us define recursively the sequence 

C Y k ’ ij ,n )k >0 as follows: for k = 0 and any (i,j) £ A\ x A 2 , we set Y°’ l3 ’ n := —Y. For k > 1, we define 
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(yk,ij,n^ z k,ij,n fU k,i3,n^ & S 2 x U 2 x H 2 (N) as the solution of the following system of BSDEs: V(i,j) G 
Ai x A 2 , 


-dY s k ’ ij ’ n = 


+n(Y s fe, b>™ — max(Y 6 
leA 1 


k-i,ij,n _ ^( S) x*>*)))-ds - Z^< n dB a - f E U*’ ij ’ n (e)N(ds, de), s<T ; 


-Trk.ij.n 

y t 


h ij (x!^ x ). 


(4.11) 

The solution of (4.11) exists and is unique since it is a decoupled multi-dimensional standard BSDE with a 
Lipschiz coefficient, noting that (YY~ 1 ’ pq ’ n )( Ptq ) & A 1 xA 2 is already given. Since n is fixed and the coefficient 
<j?i' n defined by: 


<f> ij ’ n (s, w, (/ s )( P , 9)e ^x^, Z ij , u«) := f ij (s, X*’», (y pq )( P , q )e A ^A 2 ,z ij ,u^)+n(y^- max(^-^( S , X^))Y 


is Lipscliitz w.r.t. ((y pq )(p^^A 1 xA 2 , z 12 the sequence (Y k ’ z ^ n ') k >o converges in 5 2 to as k — > 00 , 

for any i,j and n. Finally by comparison and an induction argument, used twice (with n and then with k), 
we obtain: 

yfc.y.n < y-fc,ij,n+1 j y-fc,ij,n y- 


Note that for the second inequality, we take into account of the fact that n(Y s — max(Y s — g (s, X* ,x ))) = 0 

i&Aj — u 

since g a > 0. Take now the limit w.r.t. k in the previous inequalities to obtain: 

Y ij ’ n ’ 0 < Y ij ' n+1)0 and Y ij ’ n ’° < Y. (4.12) 


Therefore there exists a progressively measurable process Y lJ, ° such that 

Y ij ’ n ’° /> Y ij ’° and Y ijfi < Y. 


Now using the monotonic limit theorem by E.Essaky ([14], Theorem 3.1) there exist (Z l i’°, f7 y, °) G 'R 2 x 
H 2 (N) and K lj{i G S 2 non-decreasing such that: 

(a) Y l k 0 belongs to S 2 and for any stopping time r, lim /• Y^ ,n ’° = Y*- 7,0 . 

n—>oo 

(b) K l k° is predictable ROLL non-decreasing, Kq’ = 0 and for any stopping time r, the sequence 
(K l J ’ n '°) n > 0 converge weakly in L 2 (P) to K Y’° ; 

(c) For any p G [1,2), 


lim E[ [ \Zi j ’ n ’° - Z ij ’°\ p ds] = 0, lim E[ [ [ \U l s j ’ n ’° - U ijfi \%n(de)ds\ = 0; 

n-s-oo J Q n-Hx J Q J E 

(d) Moreover for any (i,j) G T and s <T we have: 


Yp° = hV(XF) + /; /«(r, X*-*, (Y r w -°) (M )er, Z«>°, U? fi )dr + A^’° - 


- jf Z«>°dB r - £ / B u^°(e)n(de)dr ; 


Yb>°> ma x{Y s fc i-°- 5 (s, **■*)}; 

fee A) ~ lK 

fo T tfi' 0 ~ rnax{Y^’° - g. fc (a, X^dK^ = 0. 


(4.13) 
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The remaining of the proof is the same as the one of Theorem 3.2 in [20], pp.1623, i.e., to show that the pre¬ 
dictable process K 1 ^ 0 is continuous thanks to the non free loop property (A4). Thus (Y*- 7 ’ 0 , Z l ^°, (7*- 7,0 , RW’ 0 ) 
is a solution of (4.6) with m = 0. 

Uniqueness of the solution of (4.6) is obtained in the same way as in ([18], pp.193) or ([6], pp.122) in making 
the connection between the solutions of systems of type (4.6) and the value function of the related optimal 
switching problem. This is possible since the hypotheses on the data allow for comparison in this framework 
of Brownian-Poisson noise type (especially (Al)). 

Finally the last property of convergence stems from the following facts: (i) Y*- 7 ’’”’ 0 /*„ Y*- 7 ’°; (ii) Yj 7 ’”’ /* n 
Yj 7 ’ ; (iii) A weak version of Dini’s theorem for RCLL process (see [8], pp.202). Note that property (ii) is a 
consequence of continuity of A*- 7,0 which implies that the predictable projection of Y*- 7,0 is nothing but Yj 7 ’ 
and the same holds for Y y ’"’ 0 . 

(ii) By (4.3), (4.4) and (4.12), we obtain that the sequence of functions (v l ^ n,0 ) n >o is convergent for any 
( i,j ) € T. So let us set u l ^°{t,x) := lim„ /• v l i’ n '°{t,x). Therefore by (4.3) and (4.8), the relation (4.9) 
holds true. 

Next as previously mentionned, we can obtain the same results for arbitrary m and not only for m = 0. 
Therefore we define u l ^ m {t,x) := lim n v^' n ’ m {t,x). Those functions verify (4.9) and u 1 ^ ,m > u y,m+1 by 
(4.4). Next (« ,J,m )( i j) e ^ i x^ 2 is a viscosity solution of (4.10) since is solution of (4.1) and by 

arguing as in ([20], Theorem 4.1). Finally uniqueness in the class IT, and continuity holds true by Remark 
3.1. □ 

We now consider the increasing approximating scheme: V(i, j) € T, 

f (Y? j ’ n ,Z ij ' n ,U ij ’ n ,K ij ' n ) € S 2 x H 2 x H 2 (N) x A 2 ; 


rr = hV(x % x )+/; f^x**, {Y^ n )^ l)eT ,^ n ,ir^ n )dr - z?> n dB r 

~Sj f E Wr J ’ n (e)N(dr,de) (A^> - £«•"), a < T; 


< min{Y^" + g jt (s, X^)}, s < T, 


i n T (rr - mm {Y k r+9 jl (8’X t r})dK!r = o, 

ieAj 

where for any (i,j) G T, n > 0 and (s,y,z,u), 

l ij ’ n (s, x *■*, y, z, u) :=g ij ’-' n (s, X**, (y kl ) {k , , )eA i xA ,, *, f E u(e) e)n(de)) 


(4.14) 


xr, (y k %, l)e r,z, f E u(e)r(Xr, e)n(de)) 


+ n(y‘ J - ma x{y k ° -g (s,Xs X ))Y 
keA} ~ lK 
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The existence of (Y %3 ’ n , Z l3,n , U l3 ’ n , K l3,n )u j\ c 41 v a? is obtained thanks to Proposition 4.2 in considering 
the system of reflected BSDEs with interconnected lower obstacles associated with the data 

«))(*,which has a unique solution 
(YS i ' n ,Z i ii ’ n ,U 1 ij ' n ,K 1 i *' n ) (iJ)f : A i x a 2 and then it is enough to set ,Z* 3 ’ n ,U i3 ’ n ,K l3 ’ n ) {iJ)eAlxA 2 = 

(—hi*- 5 ’", — Zfi 3 ’ n , —U\ 3 ' n , Ki j ’ n )^ i j' )eA i xA 2 . The following is the analogous of Proposition 4.2. 

Proposition 4.3. i) For any fixed ( i,j ) £ T and n > 0 we have: 

lim E[ sup |yb>,m _ yb>| 2 ] ^ 0 and P ~ a.s., Y ij ’ n < Y ij ’ n+1 . (4.15) 

m->oo t <s<T 

(ii) There exists a unique mi x m 2 -uplet of deterministic continuous functions ( u kl ’ n )(k,i)eA 1 xA 2 in n g such 
that, for every (t,x) £ [0,T] x R fc , 

Y_f n = u ij ’ n (s,X t s ’ x ),Vs £ [t,T\. (4.16) 


Moreover, for any ( i,j) £ T and (t,x) £ [0,T] x R fc , vf 3,n (t,x) < vf 3,n+1 (t,x). Finally, (u l3,n )(ij)£A 1 X A 2 
is the unique viscosity solution in the class n g of the following system of variational inequalities with inter¬ 
connected upper obstacles. V(i,j) £ T, 


rna x{u l3,n (t,x) — min (u d,n (t,x) + gMt,x))\ —dty? 3 ’ n (t,x) — Cu l3,n (t,x) 

-g ij -’ n (t,x,(u kl ’ n (t,x))^ t l) er ,a(t,x) r D x u ij ’ n (t,x),I i j(t,x,u ii ’ n ))} =0 ; 
u ij ’ n (T,x) = h ij { x ). 

Now let us define u*- 7 ’, vf 3 , (i,j ) £ P, by: 


(4.17) 


u l3 (t,x) := lim u Z3,m (t,x), vf 3 (t,x) 

m —>00 


lim yf 3 ’ n (t,x), (t,x) £ [0,T] x R fc . 


We then have: 


Corollary 4.1. For any ( i,j ) £ A 1 x A 2 , the function id- 7 (resp. yf 3 ) is use (resp. Isc). Moreover, u 13 and 
u zj belong to II g and for any (t,x) £ [0,T] x R fe , 

uf 3 (t, x) < u lj (t, x). 

Proof. For any (- i,j ) £ A\ x A 2 , the function id 3 (resp. uf 3 ) is obtained as a decreasing (resp. increasing) 
limit of continuous functions. Therefore, it is use (resp. Isc). Next, for any (i,j) and n, m, 

u ij ’ n ’ m (t, x) < u ij ’ n ’°(t, x ), {t, x) £ [0, T\ x R fc , 

as the sequence (u l3,n ’ m ) m >o is decreasing. Thus, taking the limit as m —> 00 we obtain, 

u ij ’ n < u ij ’ n ’°. 


Now using (4.3) and (4.8), it follows that, for any t < T and s £ [ t,T ], Yf 3,n, ° = u Z3,n ’°(s,Xl’ x ) and the 
processes Y l3,n ’° converges in S 2 , as n —> 00 , to Y 13 ’ 0 which is solution of (4.6) with m = 0. Furthermore, by 
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(4.9), there exists a deterministic continuous function u lJ ’° with polynomial growth such that for any t < T 
and s G [ t,T ], Y^’ 0 = u u, 0 (s, Xl’ x ). Then taking s = t and the limit as n —» oo to obtain 

uf-’ft, x) := lim t/ J ’”(i, x) < lim u l -’’ n ’°(t,x) = u l ^’°{t, x), V(t,x) G [0, T] x R fe . 

n—> oo n—> oo 

But w *- 7,0 and u*- 7 ’" belong to II g and u*- 7 ’" < Thus u *- 7 G n^, for any ( i,j) G A\ x A 2 . The last 

inequality follows from (4.4) and the definitions of u^ and w *- 7 . In the same way one can show that ■u *- 7 G n g , 
for any ( i,j ) G Ai x A 2 . □ 

We now focus on the proof of existence of a solution for system (2.1). 

Proposition 4.4. The family (u lJ )(ij)^A 1 xA 2 is a viscosity subsolution of the system (2.1). 

Proof. First recall that for any (i, j ) G T, u *- 7 = lim m \ u^ ,m , so that ■u *- 7 is use. Moreover since (T, x) = 
h l i (x) then u P ( T , x) = (x) , Vx G R fc . 

Now let ( i,j ) G T and ( t,x ) G (0,T) x R fc be fixed. We suppose that there exists eo > 0 s.t. 

u % i(t,x) > L^[v\(t, x) + e 0 , (4-18) 

otherwise the subsolution property holds. Thanks to the decreasing convergence of (w y ’ m ) m >o to rt*- 7 , there 
exists mo such that for any m > mo, we have 

U ij ’ m (t,x) > L^[(u^ m )( P , q )^xA 2 Kt,x) + |. (4.19) 

As for any m > 0, u l ^ m and L l i\(u pq ’ m )( Ptq ) G A 1 xA 2 \ are continuous, then there exists a neighborhood 0 m of 
( t , x) such that 

U ij ’ m (t',x') > L^[(u pq m P , q )^xA 2 ](t',x’) + V{t',x') G 0 m . (4.20) 

Let now 0 be a C 1,2 -function such that < p(t,x ) = u lJ (t,x) and ■u 1 - 7 — (/> has a global strict maximum in 
(t,a:). Next let <5 > 0 and for m > 0 let ( t m ,x m ) be the global maximum of — cj) on [0,T] x B' (x,25Kp) 
(K f 3 is a bound for /3 and R'(a;, 25Kp) is the closure of 13(;r, 2 5Kp)), which exists since the function u *- 7 — (f) 
is use. But there exists a subsequence {m^} such that 

C tm k ,Xm k ) -tfc (<,*) and u l0 ’ mk {t mk ,Xm k ) -tfc U lJ (t,x). (4.21) 

Actually by Lemma 6.1 in [7], there exist a subsequence {m*,} and a sequence (t' , x' mk )k such that 

(*.*) and u l3 ’ mk (t' mk ,x' mk ) -A u l °(t,x). 

Next let us consider a convergent subsequent of ( t mk , a; mfe ), which we still denote by ( t mk , x mk ), and let (t , x) 
be its limit. Then for some ko and for k > ko we have 

u lJ (t, x) - ( j)(t,x) > limsup fc (u I - 7,mfe (t mk , x mk ) - (f(tm k ,x mk )) > lim inf k (u*- 7 > mfc (t mk , x mk ) - <f{t mk ,x mk )) 

> lim inf fe (w *- 7 ’ 7 "* 5 (C fc , J “ Wm k ,x'm k )) = u i3 (t,x) - <f(t,x) > u lj (i,x) - <j){i,x). 
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It implies that u lJ (t,x) — 4>(t,x) = u l ^(i,x) — 4>{t,x) then (t,x) = (t,x) since the maximum is strict. On the 
other hand we obviously have u l ^ rnk {t mk ,x rnk ) u l ^{t,x). Finally since this is valid for any subsequence 
of (t mk ,x mk ), then the claim follows. 

But from the subsequence {to*,} one can substract a subsequence which we still denote by {m,k} such 
that (t rnk ,x mk ) belongs to O mfc . Indeed if this is not possible one can find a subsequence {m p } of {m^} such 
that for p > 0, (t mp ,x mp ) does not belong to 0 mp , i.e., 





Then in taking the limit w.r.t. p we obtain 


where (.)* stands for the upper semi-continuous envelope. But 


(■ L lJ [{u pq ) MeA i xA 2 ])* = nra x(u kj - g )* = ma x(u kj -g ). 

keA\ — lK keA\ — lK 

Therefore we have 

u l3 (t,x) < max(u fej - 9 ik )(t,x) + ^ = L l3 [u](t,x) + ^ 
keA\ ~ lK 2 2 

which is contradictory with (4.18). Hereafter we consider this subsequence {rrik}- 

Now for k large enough: (i) (t mic , x mk ) G (0, T ) x B(x, 2 Kp5) and is the global maximum of u l3,mk — <j> in 

(0 ,T) x B(x mk ,Kp5) ; (ii) u ZJ ' mk (tm k ,x mk ) > L tJ [(u pq ' mk )( Piq)eA i xA 2 ](t mk ,x mk ). As u lJ ’ mk is a subsolution 

of (4.10), then by Proposition 5.1 - Remark 5.1 in Appendix, we have 

dtfiitmk > Xmk ) — £-4 > {tmk > X mk ) Ijj (fmk ) X mk > < / > ) > Xmk > D x (j)(t mk , X mk \u 3 k ) 

+m k (u x k (t mk , x mk ) — ' k {t mk ,x mk ) + gji(t mk ,x rnk ))) + (4 22) 

— 1? i [^mfc > Xm k i ’ k (t‘m k i Xm k ) ) (p,l)G A 1 xA 2 : ^{^m k > Xm k ) B) x (j)(t rnk > Xm k ) > 
lij (tm k ,X rnk ,<t>) + IA (t mk ,X rnk ,U , X rnk )}. 

From which we deduce, in dividing both hand-sides of (4.22) by mk and then taking the limit as k -A oo, 
that 

Cfc — (w X k (tm k > Xm k ) min (u k (t-mk > Xm k ) 4“ gjl (tm k : X mk ))) ^fc 0. 

Next fix k 0 and let k > ko. As the sequence {u l3,m ) m is decreasing then 

k {t mk ,x mk ) < nrin ;g ^j (u ’ k (t mk ,x mk ) gji{t m ki x mk )) + e k 

— (^ ’ k ° (tm k > 2-mfc ) + djli^mk i X mk )) + €fc 

Take the limit w.r.t fc, using continuity of u ll,mk o then send kg to +oo to obtain: 


u l3 {t,x) < min (u ll (t,x) + gji(t,x)). 
leA{ 


Next there exists a subsequence of {to*,} (which we still denote by {rtik}) such that: 
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(i) V(p,Z) e A* x Afj, (■ u pl ’ mk (t mk ,x mk ))k is convergent and then \im k u pl ’ mk (t mk ,x mk ) < u pl (t,x) ; 
(u) Ig {tm k 1 x m k ) 4 >) ->fc Ig(t,x,<j>) and llf(tm k ,x mk , </>)) ->k i]f(t,x,4>)i 

(iii) By Fatou’s Lemma, limsup k I$(t mk ,x mk , D x (j>(t mk ,x mk ),u i:l ’ mk ) < lf{t,x,D x <j){t,x),u i:l ) and 
limsup k lff(tm k ,X mk ,U ij ’ mk ))) < iff (t, X, u ij ). 

Let us now set: 

•— 9 j 5 ^Tn k 5 ’ k (trrik 5 Xmy. ))(p,/)G-Ai X A 2 i ^(j'rrik 5 ) Dxtfi(j'm k •> %m k )? 

(*m fe , *m fc , 0) + iff (tm k , ^m fc , U ij ’ mk )] 

9 ^ \pmk 5 ^rrifc 5 ’ fc (^m/.; ))(pi)GA 1 X A 2 ) 5 %/-) Dx^i^mk 5 x m k )5 lij 0 ) Hb L—’ 

Then, by linearizing w.r.t. <?, there exists a non-negative bounded quantity S 2 such that 


Afc — ^2 X (Xj {tm k ) x m k ) </>) + Zj (in 


iff 


< Cij x (Jx {t mk ,x mk ,(j)) + 1^ (t 


, jj'-)'™*) — ijf(t, x, 4>) 


m k 5 ^mfc 


lff{t,x,u 13 ) 
lff(t,x,u lJ )) + 


where Cij is the Lipschitz constant of . Therefore, with (ii)-(iii) above, we have that lim sup /c A* < 0. 
Going back now to (4.22), and take the limit superior w.r.t. k to get: 

-d t <j>(t,x) -£<j>{t,x) < Ig(t, x, (j>) + lf(t,x, D x (/>(t,x),u i3 )+ 

9 ij [t, x , (■ u pl {t , x)) MeT ,a{t, x) T D x <j){t, x),lff(t, x, <j>) + lff(t , a;, u ij )\. 

But u^(t,x) — (f>(t,x) and uP < 0, then lf(t,x, D x (f>(t,x),u l:i ) < lf(t, x, D x cj)(t, x), <j>) and iff {t,x,u'i) < 
iff it, x, 4>) ■ Plugging now this inequality in the previous one to obtain 


-d t 4>(t,x) -£<j)(t,x) -I{t,x,(j>) - g l: >[t,x,(u pl (t,x)) MeT ,(j(t,x) T D x u i: >{t,x),Iij(t,x,(j))\ < 0 . 
Therefore is a viscosity subsolution of 

min{(w - L^[(u kl ) {kil)er ])(t,x);max{(w - U ij [{u kl ) {ktl)er ])(t,x); 

-d t w{t,x) - Cw(t,x) - g lJ (t,x, [(u pl (t,x))^ l ) er ^i, j ),w],a{t,x) T D x w{t,x),Iij{t,x,w))}} = 0 ; 
w(T, x) = h l i{x). 

As (■ i,j ) in T is arbitrary then {u z: ’)(ij)^A 1 xA 2 is a viscosity subsolution of (2.1). □ 

Proposition 4.5. Let mo be fixed in N. Then the family (u l3,rn °)^i t j) e A 1 xA 2 is a viscosity supersolution of 
system (2.1). 

Proof. Recall that (? ,J ’ m °, U l3 ’ m °, K l I m °)^ j)^a 1 xA 2 solves the system of reflected BSDEs (4.6). 

Therefore if we set K l J' mo f“ := ?n 0 J^(Yf 3 ’ mo — min(l^ i,m ° + gji(r, Xf x ))) + dr, s <T, then 

leA2 i 

(yii.mo ) ) ^ij>„ ~)(i,j)eA'xA 2 is a solution of the following system of reflected BSDEs 
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with bilateral interconnected obstacles: for any (i,j) £ F and s <T, 


yb.mo =h i J(X^ x ) + J s fi (r, X*’ x ,(Y r kl ’ m °) {ktl)er , Z?’™ 0 ,U?. j ’ mo )dr - / s Z l r 3 ’ m °dB r 

- if Ie U^ ,mo (e)N(dr, de) + Kf m ° - K 

max{? s fe ^ m ° - 3., (s, X*-*)} < Y^’ mo < V m i n (V?h m ° + m.,(s, W’ x )); 

fceAj — lk IgA 2 j 

fo(Xs j ’ m ° - -5. fc ( S ,X*> x )})d^™o - 


fo( Y s ij,mo -Yji’ mo \/m.m(Y s il ' mo +g jl {s,X t s ’ x )))dK i 7 ?’ m ° ~ = 0 


(4.23) 

On the other hand we know by (4.9) that there exist deterministic continuous functions (u lJ ’ m °)rij)^ A i xA 2 
in n 9 such that, for every t <T, 


Y s ij ' m 0 = u ij ’ mo (s,Xl’ x ), s £ [t, T}. 


Then using a result by Harraj et al. [21] we deduce that for any (i, j) £ F, u l3 ' m ° is a viscosity solution of 
the following IPDE: 

min {u ij ’ mo (t,x) - max fcgA i {u kj ’ m °(t,x) ~ 9 ik (t,x)}; 

ma x{u^ ,m ° ( t, x) — « y,m0 ( t , x) V min /gA 2 ( u ll ’ m ° ( t , x) + gji ( t , a;)); —dtu' l3 ’ mo ( t, x) — Cu l3,rn ° ( t , a;) 

(f. (u kl ’ m ° (t, a;)) (/M)gA 1 x A 2 , cr(t, x) T D x u ij ’ m ° (t, x),l l3 (t, x, u ij ’ m °))}} = 0; 
u ij ’ mo (T,x) = h ij (x). 

But 


u^’ mo (t,x) — u l 3 ’ m °(t, x) V mini &A 2(u il,rn ° (t, x) +9ji(t,x)) < u l 3 ' m °(t,x) — mini GA 2(u zl ’ m ° (t, x) +gji(t,x)) 
Therefore u l3 ' m 0 is a supersolution of 

min{u^> mo (t,x) - max keA i{u kj ’ mo (t,x) ~ 9 ik (t,x)}\ 

ma x{u l 3 ’ m °(t,x) — niin igA 2 (u ll,m °(t,x) + 9 ji{t,x)); — d t u l 3 ’ m °{t, a:) — £u' l 3 ’ m °(t, x) 

- 9 ij (t, x, (u kl ' m ° ( t, x)) {ktl]eA i xA 2, a(t, a :) T D x u ij ’ m ° (t, x), % (t, x, #"'«))}} = 0 ; 
u ij ’ m °(T,x) = h ij {x). 

As (i,j) is arbiratry then (« I -’’ m °)(ij)gA 1 xA 2 is a viscosity supersolution of system (2.1). □ 

Consider now the set U mo defined as follows: 

U 7 n 0 = {u := (« y )(ij) € A x xA 2 s.t. u is a subsolution of ( 2 . 1 ) and V(i,j) £ A 1 x A 2 ,u 13 < u 13 < u l3 ’ mo }. 
W mo is not empty since it contains (« y )(ij)eAi X A 2 - Next for (t,x) £ [0,T] x R fc and ( i,j ) £ T, let us set: 

mo u l 3 {t,x) = sup{u y (t,x), (u kl ) {k ' l)eA i xA 2 £U mo }. 


We are now ready to give the main result of this paper: 
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Theorem 4.1. The family ( rn °u l3 )uj) &A i xA 2 does not depend on mo and is the unique continuous viscosity 
solution in the class II g of the system (2.1). Moreover m °u 13 = u 13 for any ( i,j ) G A 1 x A 2 . 

Proof. Firs note that w.l.o.g we assume that for any (i,j) G A 1 x A 2 , the function 

l/eR4 gij(t,x, > y]i z i o) is also non-decreasing when the other variables are fixed. 

To begin with, note that for any ( i,j ) G A 1 x A 2 , u 13 < TOo u u < u l3,rn °. Since u 13 and u 4 - 7 ’ 4 " 0 are of 
polynomial growth, then { m °u l 3 )^j\ eA i xA2 belongs also to II g . The remaining of the proof is divided into 
two steps and for ease of notation, we denote simply by (u 4 - 7 ^j,p e p as no confusion is possible. 


Step 1: First we show that (u 13 )^j') eA i xA 2 is a subsolution of (2.1). As u 13 < it 4 - 7 < u Z3 ’ m ° then 
w 4 - 7 < u 13 ’* < u y,m ° since tt 4 - 7 is use and is continuous. Therefore, for any x G R fe , since u l3 (T,x) = 

u l3 ' m °(T,x) = h Z3 (x), we have u Z3 ’*(T,x) = h Z3 (x). 

Next let {u l - i )(ij)eA 1 xA 2 be an arbitrary element of U mo and let (i,j) be fixed. Let (t, x) G (0,T) x R fc 
and <j) G C 1,2 ([0,T] x R fc ) such that u l3 '*(t,x) = 4>{t,x) and u 13 '* < (j>. Then 

min{(u y ’* - L ZJ [{u kl ’*)(k,l)eA 1 xA 2 ]){t, x)\ max-J^w 4 - 7 ’* - U l3 [(u kl ’*)( k> i )eA i xA 2 ])(t,x)i 

-d t (t>(t, x) - C<j>{t , x) - g ij {t, x, ( u kl ’*(t , x))( ki i) eA i xA 2 , a(t, x) T D x <j)(t, x), I zj (t, x, 0 ))}} < 0. 

By definition, for any (fc,Z) G T, u“ < and then u W ’* < Using now the monotonicity property 

(A2), we obtain 

min{(u^'>* - L i3 [(u kl ’*)( ky i )eA i xA2 ])(t,x);max{(u ij ^ _ U ij [(u kl ’*) {kt i )eA i xA 2 ])(t,x); 

-d t 4>(t,x) -C<j>(t,x) - g i3 {t,x,[(u kk *) {kl) ^ T - {i ,o,u i3 ’%t,x),a(t,x) T D x cj)(t,x),lij(t,x,(t>))}) < 0. 

This means that u 13 is a subsolution of the following equation: 

f min{(w - L l J[{u kl ’*)( k ,i )&A i xA 2 })(t,x);ma,x{(w - U ij [(u kl '*) (kt i )eA i xA 2 ])(t, x); 

< —d t w(t, x) - Cw(t, x) - g ZJ (t, x, [(« fci ’*) (fc ,; )er - ( M), io](t, x), aft, x) r D x w(t, x), x, to))}} = 0; 

I w(T,x) = h Z3 (x). 

Consequently, by a result by Barles-Imbert ([4], Theorem 2, pp.577), w 4 - 7 j s a subsolution of 

f min{(w - L l J[(u kl ’*) {k ' l)eA i xA 2 ])(t,x);max{(w - U ij [(u kl ’*)^ l)eA i xA 2 ])(t,x); 

< —d t w(t, x) - Cw(t, x) - g ZJ (t, x, [(« fci, *) (fc ,i) er -(M), x), cr(t, x) T D x w[t, x), /^(t, x,«;))}} = 0; 

I w(T, x) = h Z3 (x). 

As (i, j) in T is arbitrary then (u Z3 )^jy €A i xA 2 is a subsolution of (2.1). 

Step 2: We will now show, by contradiction, that (u Z3 )uj\ eA i xA 2 is a supersolution of (2.1). First note 
that for any (i,j) G A 1 x A 2 , u) 3 = yfj < -rfi < u 4 - 7 < jT 3,m ° = u Z3 ’ m °, since u 4 - 7 ’ 4740 is continuous and uf 3 j s 
Isc. Therefore, for any x G R fc , since u Z3 (T,x) = h Z3 (x) = U 47 ’™ 0 (T,x), it holds, uJ(T,x) = h l3 (x). 
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The rest of the proof is rather classical and can be read e.g. in [5] up to some adaptations. However we 
give it for completeness. So suppose that there exist (i,j) £ A 1 x A 2 , (t, x) £ (0, T) x H k and a C 1 —function 
<j) such that uj(t,x) = cp(t,x), u 2 (s,y) > <p(s,y) in (0,T) x R fc — {(i,a:)} and 

min{((p - L ij [(u kl )( k: l )eA i xA 2 ])(t,x)-,max{((j) - U ij [(u kl ) {kt i )eA i xA2 ]){t,x); 

-dt<t>(t, x) - C(j>(t, x) - g^{t, x, [(u kl ) (fc , i)er -u,n, <P](t, x),a{t, x) T D x <p(t, x), I^t, x, </>))}} < 0. 


Then by continuity of the equation, continuity and monotonicity of gA and lower semi-continuity of u kl , 
there exist two constants > 0 and <!>i > 0 such that: V(s, y) £ B((t, x), S i) and 0 < e < q we have: 

min{(^ e - L^[(u kl ) {ktl)eA i xA 2 ])(s,y);max{((f> e - U ij [(u kl )( k>l)eA i xA 2 ])(s, y)\ 

-dt^eis, y) - C<j> e (s, y) - g tJ (s, y, [(«*0(M)er-(M), </> e ](s, y), a(s, y) T D x <p e (s, y), y, (p e ))}} < 0 

(4.24) 

where <p e = <p + e. Next since (t, x) is a strict minimum of up — <p , there are constants 0 < ti and 0 < 62 < <5i 
such that up — <p > e 2 on 0B((t, x), 62 )- Now let us set 63 = min(ei, e 2 ) and let us define {w kl )( k n eA i xA 2 as 
follows: 


w kl = u ki ,* if ^ ^ (i } j) and w ij 


max(<j> + € 3 , it®- 7 ’*) on B((t,x),S 2 ) C (0, T) x R fc ; 
u l i’* elsewhere . 


Then (w kl )^ eA % xA 2 belongs to n 9 and is a viscosity subsolution of (2.1). Indeed, first note that for any 
(k,l) £ T, w kl is use and w kl (T,x) = u kl, *(T, x) = h kl (x). Next let (s,y) £ (0 ,T) x R fc . If (s,y) does not 
belong to B((t,x),S 2 ) then the subsolution property stems from the one of {u kl ’*)( k n eA i xA 2 . Assume now 
that (s,y) £ B((t,x), 52 )- If (k,l) ^ (i,j), then the subsolution property stems from the one of u kl ’*, in 
taking into account of re ®- 7 > it®- 7 ’*, w kltl = u kl11 ’* if 7 ^ (i,j) and monotonicity of g kl . Finally let us 

deal with the case when ( k,l ) = (i, j). Let ^ be a C 1,2 -function such that w l ^{s,y) = i(>{s,y) and xp — tu ®- 7 
has a strict global minimum in (s,y) £ (0, T) x R fe . If w l i{s,y) = u l i'*(s,y) then 


min{(V> - L ij [(u kl ’*)^ k! i) eA i xA 2 ])(s,y);max{(ip - U ij [(u kl ’*)( k j )eA i xA 2 ])(s,y); 

-d t xp(t,x) - Cp’{s,y) - g lJ (s,y, [(u kl ’*) {k j )er -a,j), xp](s,y),a(s,y) T D x xp(s,y),Iij(s,y, xp))}} < 0 

since w®- 7 > u lJ '* and then xp — it*- 7 ’* has a strict global minimum in ( s,y ), (ip — u lJ ’*)(s,y) = 0 and 
{ ukl ’*)(k,l)eA 1 xA 2 is a subsolution of system (2.1). Thus, 

min{(V’ - L^[(w kl ’*) {kt i )eA i xA 2 ])(s,y);max{(xp - U ij [(w kl ’*)( kt Q eA i xA ?])(s,y); 

-d t ip{t, x) - Cxp(s, y) - g 13 (s, y, [(w fci ’*)( M ) er -o .:0 , ip] (s, y), er(s, y) T D x xp{s, y), (s, y, V 7 ))}} < 0 . 

(4.25) 

Finally if w lJ (s, y) = <p{s, y) + 63 then <p(s, y) + e 3 = xp(s, y) and <p + €3 < xp on B((t, x),$ 2 )- It implies that 
d t (j>(s, y) = d t xp{s,y),D x <p(s,y) = D x xp{s,y) and D 2 xx <p(s,y) < D 2 xx xp(s,y) 


and by (4.24) we deduce that tu®- 7 verifies (4.25). Therefore in®- 7 satisfies the subsolution property and 
( 7 a fei )(fc,z)e- 4 1 xA 2 is a viscosity subsolution of (2.1). But we have, 

x) > ma x{<p(t, x) + 63 , u z J(t, x)} = (pit, x) + €3 = u l J(t, x) + € 3 . 
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Thus there exists (to,xo) G (0 ,T) x R fc such that w l 3 (to,Xo) > u t 3 {ta,xo). But this is contradictory to the 
definition of u 13 . Therefore (« w )(M)eA 1 xA 2 is a supersolution of (2.1) and the proof is complete. 

Now, by Corollary (3.1), { m °u l 3 )/ i j\ eA i xA 2 (he. ( ut; ')(ij)eA 1 xA 2 ) does not depend on ?n 0 since the 
solution of (2.1) is the unique. On the other hand for any ( i,j ), we have 

u ij < u ij := m °u ij < u ij ’ m ° 

and in taking the limit wrt mo we obtain u l ° = u 13 , for any (i,j) G A\ x A 2 . □ 

As a by-product of the above construction we have the following result related to the limit of the increasing 
scheme: 

Theorem 4.2. The family (u l 3 )^j-\ eA i xA 2 continuous and of polynomial growth and is the unique viscosity 
solution in n g of the max-min problem, i.e., for any ( i,j ) G T, 

{ max{(■ v 13 — U 13 [€>]) (t, x) ; min{ ( v 13 — L 13 [w]) (t, x); 

-d t v ij (t,x ) - Cv ij (t,x) x, (v kl (t,x))( k g )eA i xA 2,cr{t,x) T D x v i 3 (t,x),r i 3 (t,x,v i3 ))}} = 0 ; 

v ij (T,x) = h i 3 (x). 

(4.26) 

Proof. Actually in considering the opposite of the increasing scheme defined in (4.14), which becomes a 
decreasing one, we obtain that (—u l 3 )(ij)£j ii x a 2 is continuous and of polynomial growth and is the unique 
viscosity solution in II g of the following system: V(i, j) G A 1 x A 2 , 

{ min{(yb - max keA 2{v ik - g k j})(t,x);max{(v ij -min ieA i{v y +g. l })(t,x); 

-d t v ij (t,x) - £v ij (t,x) +g i 3 (t,x, (-v kl (t,x)) (ki i )eA i xA 2,a(t,x) T D x (-v i 3 )(t,x),-I l:) (t,x,v i3 ))}} = 0 ; 
v iJ (T,x) = -h ij (x). 

Using now a result by G.Barles ([2], pp.18) we obtain that {u i: ’)(i,j)eA 1 xA 2 is the unique viscosity solution 
in n g of (4.26). □ 
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5 Appendix : Alternative definition of the viscosity solution of 
system (2.1) 

The following result inspired by the work by Barles-Imbert [4] is another definition of the viscosity solution 
of system (2.1). We do not give its proof since it is an adaptation of the one given in ([20], Proposition 5.2, 
pp.1656) as the function q i —> g lJ (t, x, y, z, q) is non-decreasing, /? is a bounded function, 7 y is non-negative 
which then imply Iij(t,x,<f>) < Iij(t,x,ip), iff (t, x, (j>) < lff(t,x,ip) and lff(t,x,<j>) < iff (t,x,tf>) for any 
<j) < if such that <j>(t,x) = ip(t,x) = u^(t, x) (<5 > 0 and (i,j) G T are fixed). 

Proposition 5.1. A function u = ®))(i,j)er : [0,T] x R fe —> R, m i xm 2 such that for any ( i,j) G P, 

u l i G n g is Isc (resp. use) is a viscosity supersolution (resp. subsolution) of (2.1) if: 

(i) u* J (T, xo) > (resp. <) h v (x o), Vxq G R fe ; 

(ii) For any (to,Xo) G (0,T) x R fc , 6 G (0,1) and a function <f G C 1,2 ([0,T] x R fe ) such that u l i(t 0 ,xo) = 
4>(to,xo) and u — 4> has a global minimum (resp. maximum) at ( to,xo ) on (0,T) x B(xo,6Kp) where Kg 
is the bound of (3 (see the first inequality of (AO)-(ii)), we have: 

{ min{ (u 13 - L l:1 [u])(t 0 ,x 0 );max{(«' 3 - U^[u\)(t 0 ,x 0 )-, 

~d t <l>(to,xo) ~ b(t 0 ,x 0 ) T D x (/)(to,Xo) - \Tr(acr T (t 0 ,x 0 )Dl x (j)(t 0 ,Xo)) - Ig(t 0 ,x 0 ,<j>) - lf(t 0 ,Xo,D x (/),u i:l ) 
-g ij (to,xo,(u kl (to,x 0 ))( k j )eA ixA 2 ,(r(to,xo) r D x (l>(to,xo),ll/(to,xo,(p) + lff(to,xo,u ij ))}} > (resp. <)0. 


Remark 5.1. In taking gji = +00 (resp. g. k = + 00 ) for any j,l G A 2 (resp. i,k G A\) we obtain an 
alternative definition of the viscosity solution of the system of variational inequalities with interconnected 
lower (resp. upper) obstacles. □ 




